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We study the dynamics of three particles in a finite interval, in which two light particles are sep- 
arated by a heavy "piston" , with elastic collisions between particles but inelastic collisions between 
the light particles and the interval ends. A symmetry breaking occurs in which the piston migrates 
near one end of the interval and performs small-amplitude periodic oscillations on a logarithmic 
time scale. The properties of this dissipative limit cycle can be understood simply in terms of a 
effective restitution coefficient picture. Many dynamical features of the three-particle system closely 
resemble those of the many-body inelastic piston problem. 
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I. INTRODUCTION 

In the preceding paper, denoted as HR p| , we discussed 
the collision dynamics of an elastic three-particle system 
on a finite interval that consists of a massive particle — a 
piston — that separates two lighter particles. The mo- 
tivation for studying this idealized system was to shed 
light on the enigmatic piston problem 0, where a gas- 
filled container is divided into two compartments by a 
heavy but freely moving piston. When the gases in each 
compartment have different initial thermodynamic states 
and when the piston moves without friction, the approach 
to equilibrium is unexpectedly complex and still incom- 
pletely understood 0, 0,1101 • 

As discussed in HR, some of the rich phenomenology 
of the piston problem can be captured by the much sim- 
pler three-particle system in a finite interval. To under- 
stand the evolution of the latter system, it proved con- 
venient to map the trajectories of the 3 particles on the 
line onto an equivalent elastic billiard particle that moves 
within a highly skewed tetrahcdral region, with the spec- 
ular reflection whenever the billiard hits the tetrahedron 
boundaries 0, IS U E3> El • From this simple geometrical 
mapping, we deduced several anomalous dynamical prop- 
erties of the three-particle system, such as the power-law 
distribution of time intervals for the piston to make suc- 
cessive crossings of the interval midpoint. 

In the inelastic piston problem, the collisions between 
the constituent particles in the gas are inelastic, so that 
each gas undergoes inelastic collapse if either the number 
of particles is sufficiently large or the restitution coeffi- 
cient is sufficiently small. Recent work by Brito et al. 01 
has again discovered surprisingly rich dynamics, very dif- 
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ferent in character from the elastic case, in which one of 
the gases cools more quickly and gets compressed into a 
solid by the piston. An even stranger feature is that this 
compression is not monotonic, but rather the piston has 
superimposed oscillations whose period grows exponen- 
tially with time. Thus the cooling of the inelastic piston 
problem is much richer than that of the classical inelastic 
gas problem [l3| . 

Given the complex behavior exhibited by the many- 
body piston system, we are again led to investigate a 
simpler alternative: a three-particle system in the unit 
interval that contains a heavy piston that lies between 
two light particles. Collisions between light particles and 
the ends of the interval (henceforth termed walls) are in- 
elastic, to mimic the many-particle piston problem when 
the gases are inelastic, while the collisions between the 
particles and the piston are elastic. 

When the light particles have the same initial energy 
but nonsymmetric positions, one light particle loses en- 
ergy more quickly than the other. As a consequence, the 
piston migrates to the wall that is closer to the cooler 
light particle. Somewhat unexpectedly, a typical system 
eventually falls into a periodic state on a logarithmic 
time scale where the piston undergoes small-amplitude 
oscillations near one wall with a constant period in hit, 
while the light particles undergo complementary oscil- 
latory motions. We term this phenomenon as the log- 
periodic state. Another intriguing aspect of the three- 
particle system is that it closely mirrors the time evo- 
lution in the many-particle inelastic piston system [l2| . 
Thus we are able to understand features of the many- 
body problem in terms of simple physical pictures that 
arise from studying the three-particle system on the in- 
terval. 

In the next section, we describe the two basic dynami- 
cal features of the three-particle system, namely, the ini- 
tial symmetry breaking and the log-periodic state. We 
then give a macroscopic description of the collapse pro- 
cess and the subsequent oscillatory motion of the piston 
in Sec. III. Finally, in Sec. IV, we develop an effective 
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FIG. 1: (Color online) Piston position X2(t) versus tona log- 
arithmic time scale for rri2 = 100 and r = 0.9. The solid curve 
is the simulation result while the dashed curve is the predic- 
tion from the macroscopic equations of motion, Eqs. (0-10. 
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FIG. 2: (Color online) Magnification of the long time evolu- 
tion in Fig. Bottom: The log-periodic state. Top: Detail 
of a "rattling" collision sequence between the piston and the 
trapped light particle. 



restitution coefficient description for the particle colli- 
sions that accounts for many of our observations. Various 
calculational details are given in an appendix. 



II. BASIC PHENOMENOLOGY 
A. Symmetry breaking 

For the many-body system in which the gases on either 
side of the piston are inelastic and have identical macro- 
scopic initial conditions, Brito et al. found an insta- 
bility in which one of the gases cools more rapidly and the 
piston ultimately compresses the cooler gas into a solid. 
While such an instability seems intuitively plausible, an 
unexpected feature is that the piston moves nonmonoton- 
ically during this cooling, with regular oscillations that 
are periodic on a logarithmic time scale. In this section, 
we show that much of this phenomenology also arises in 
the idealized three-particle system on the unit interval. 

The particles are located at x\, X2, and X3, with < 
x\ <x-i<xj,<\. The light particles, with masses vn\ — 
TO3 = 1 and locations x\ and 2:3, collide elastically with a 
massive piston with mass m 2 3> 1 at X2 and inelastically 
with the walls. Thus a light particle that hits a wall with 
speed v — \J2E is reflected with speed rv, where r G [0, 1] 
is the restitution coefficient. The energy change in this 
collision is AE = -£7(1 - r 2 ) < 0. 

Figure ^ shows a representative result for the piston 
position X2 (t) versus i on a logarithmic scale for the 
case ra 2 = 100 and r = 0.9. The initial velocities are 
(vi(0),t> 2 (0),t> 3 (0)) = (1,0,-1) so that two light par- 
ticles approach the piston with equal and opposite ve- 
locities. Thus the system initially has zero momentum 
and total energy £7 = 1. The initial positions of the 



light particles were chosen uniformly in (0,1/2) and in 
(1/2,1); for the example of Fig. □ (ii(0), x 2 (0), x 3 (0)) = 
(0.083 25,0.5,0.862 83). As a result, the first collision 
is between the piston and particle 3. This small initial 
asymmetry eventually drives the piston from oscillations 
about X2 = 1/2 to the nonsymmetric long-time behavior 
depicted in Fig. It bears emphasizing that the phe- 
nomenology of the three-particle system up to approxi- 
mately 10 5 time steps is qualitatively similar to that of 
the many-particle inelastic piston problem |12J. 

In the long time limit, the piston migrates close to one 
of the walls. Which of the two walls is selected is deter- 
mined by the identity of the first collision. When the pis- 
ton is initially located at x 2 = 1 /2 and the two particles 
approach with equal and opposite velocities, the piston 
is driven to the right wall if the first collision occurs with 
its right neighbor and vice versa. The light particle that 
first hits the piston then collides earliest with the wall 
and begins cooling earlier. This fact leads to the piston 
eventually compressing the particle that experiences the 
first collision with a wall. 



B. The log-periodic state 

Numerically, we find that the three-particle system 
asymptotically falls into a log-periodic state-where the 
piston undergoes small-amplitude oscillations with a con- 
stant period in hit — for almost all initial conditions. In 
this state, one of the light particles is trapped in a small 
gap between the piston and the wall (Fig. EJ, while the 
other light particle has most of the energy and travels 
over almost the entire interval. 

During these oscillations, the light particle that is com- 
pressed by the piston performs a sequence of violent 
rattlings each time the piston approaches the wall that 
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FIG. 3: (Color online) Plot of r as a function of m2 on a 
semilogarithmic scale, for restitution coefficients r = 0.875, 
0.8875, 0.9, 0.9125, 0.925, 0.9375, and 0.95 (bottom to top). 
The inset shows the data collapse of r for different values of 
r using 7 = 2.1 and r\ — 2.6. 

eventually reflect the piston from the wall (top panel in 
Fig. 0) . The piston then collides with the other light par- 
ticle whose energy is nearly equal to that of the entire 
system and whose momentum is comparable in magni- 
tude to that of the piston. After this collision, the piston 
is reflected back toward the nearer wall and the rattling 
sequence with the trapped light particle begins anew. 

Generally this long-time state has a one-cycle periodic- 
ity in which the position of the piston recurs at each max- 
imum of its oscillation cycle (Fig. [2J> ■ However, for piston 
mass m2 less that a r-dependent threshold mass /U*(r), 
we empirically find that the asymptotic state can be a 
two-cycle, three-cycle, etc., with lower cycles more likely 
to occur than high cycles. Conversely, for 7712 greater 
than an upper threshold 

., M _ (l + r)(l+r + Vr) + 4r 

VcV ) - (1-r) 2 ' W 

inelastic collapse occurs, where the piston ultimately 
sticks to a wall (see the Appendix for the derivation of 
fi c ). For the purposes of the present discussion, we are in- 
terested in the case where 777,2 is in the range \pt( r )t A t c( r )] 
so that the system falls into a one-cycle log-periodic state. 

This state may be characterized by the relaxation time 
r(m2,r) until the piston settles into the log-periodic mo- 
tion and the amplitude, A(ni2,r), and period on a loga- 
rithmic time scale, A(m2,r), of the ensuing oscillations. 
The latter is defined via tk+i — e A tk, where tk and tk+i 
are the times for two consecutive maxima of x 2 (t) in the 
final state (bottom panel in Fig. [2J) ■ 

Figure |3 shows the relaxation time r as a function of 
m 2 for m2 in the range \jjbt , /i c ] for representative values 
of r. We expect that the oscillatory regime is reached 
more quickly for larger m 2 since energy is more quickly 
dissipated when the piston is heavier, as confirmed by 




FIG. 4: (Color online) The amplitude of the log-periodic os- 
cillations, A, as a function of the piston mass on a semiloga- 
rithmic scale for the same r values as in Fig. El (data shifting 
to the right for increasing r). The curves are the predictions 
from Eq. iH2t . based on an effective restitution coefficient pic- 
ture (see text). Inset: Data collapse of A for different values 
of r and 7 = 2.1. 

the data. We also find that r decays exponentially with 
m 2 , that is, r ~ exp[— m 2 //i(r)], with a characteristic 
mass scale fi(r) that is nearly equal to the threshold mass 
/it(r). Both /i and [it numerically scale as (1 — 7') -7 for 
r close to 1, with 7 w 2.1. This is slightly larger than 
the anticipated exponent value of 2 that is based on the 
hypothesis that there should be only one characteristic 
mass that scales as /i c ~ (1 — r)~ 2 in the limit r->l from 
Eq. QJ. We attribute the discrepancy in 7 to corrections 
to scaling; the largest restitution coefficient r = 0.95 that 
is practical to study is still not very close to 1. 

When the piston is in the log-periodic state, the am- 
plitude A is a monotonically decreasing function of 7712 
and vanishes as 777,2 - ► Mc(?"); signaling the onset of 
inelastic collapse (Fig. 0J. The logarithmic period of 
the oscillations A (not shown) scales approximately as 
A(m2,r) « (1 — r) and depends weakly on m 2 . 

Because these three characteristics of the oscillations — 
t, A, and A — seem to be governed by the same mass 
scale, we anticipate that data collapse will occur. Em- 
pirically, we find that these quantities are consistent with 
the scaling forms 

r(m 2 ,r) ~ (1 - r)' 11 $ r [m 2 (l - r) 7 ] , 
A(m 2 ,r) ~ (1-r) 7 <5> A [m 2 (l - r^} , (2) 
A(m 2 ,r) - (1 - r) $ A [m 2 (l - r) 7 ] , 

with 7 2.1. Additionally, 77 ~ 2.6 is an apparently 
independent exponent that characterizes the relaxation 
time t, while $ T , <&^, and $a are scaling functions. The 
insets to Figs. and 0] show that the data collapse for r 
and A is quite good. 

Although the asymptotic state of the system is periodic 
on a logarithmic time scale, we emphasize that the total 
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FIG. 5: (Color online) Energy of each particle as a function 
of time on a double logarithmic scale for the system depicted 
in Fig. [I] 



energy of the system, E(t), continues to dissipate due to 
inelastic collisions with the walls. At a coarse-grained 
level, we recover Haff's law [l4| E(t) ~ V 2 (upper curve 
in Fig.[SJ), as expected. Notice that for the specific exam- 
ple being studied (in which the piston compresses parti- 
cle 3), E(t) w E\(t). At a finer time scale, however, 
E(t) undergoes a sequence of steps and almost constant 
plateaus. The largest drop in energy occurs when the 
more energetic light particle collides with the wall, while 
the rattling dynamics between the piston and the other 
light particle leads to a small decrease in the energy of 
the system. 



III. MACROSCOPIC DESCRIPTION 

We can understand the initial instability of the pis- 
ton in terms of macroscopic equations of motion 1 1 l2L Il4| . 
The macroscopic approach given here ostensibly applies 
for any value of the restitution coefficient r for fixed pis- 
ton mass m2, or equivalently for any to 2 for fixed r. In 
particular, this approach correctly describes the initial 
instability of the inelastic piston for any value of the pa- 
rameters r and mi- The only feature that the macro- 
scopic approach fails to describe is the final state for pis- 
ton masses below the critical mass, mi < 

According to the macroscopic description, the energies 
of the light particles change in time by two processes. 
First, the energy decreases due to inelastic collisions with 
the walls. This cooling is macroscopically described by 
Haff's law, in which the energy change is proportional 
to the particle energy and the number of collisions per 
unit time. Thus dEi(t)\ co \\ = —Ei(t)(l — r 2 )n co n(t)dt, 
with i — 1,3. The collision rate may be approximated by 
n co \i(t) w y/2Ei{t) / ii(t) , where ^/2E i {t) is the thermal 
velocity and £i(t) is the length of the region available for 
particle i; thus £i(t) = x%{t) and ta{t) = X — x-2.it)- On 



the other hand, the energies of the light particles also 
change because of compression or expansion by the pis- 
ton. The macroscopic equation describing this process is 
dEi(t) | piston = —Pidli, i = 1,3, where Pj is the pressure 
exerted on particle i, and d£i is the length change 

Assuming the ideal gas law Pj = Ti/ii (Boltzmann's 
constant is set to 1), and writing Tj = 2Ei, we obtain 

^ = -2^-V2(l-r»)^, (3) 

dt X2 X2 

^ = 2E 3 -^-V2(l-r>)^. (4) 

dt 1 — X2 1 — X2 

This is essentially the approach of Haff [T4| for the inelas- 
tic gas, and it was also adopted for the inelastic many- 
particle piston problem of Brito et al. 01 • The force 
exerted on the piston is given by the pressure difference, 
P\ — P2 , so that the macroscopic equation of motion for 
the piston is 

d 2 x 2 2E t 2E 3 

m2 ^2~ = 1 ■ ( 5 

dt 1 x 2 1 — x 2 

Equations describe the evolution of the three- 

particle system on a coarse-grained time scale and they 
are the analogs of the equations derived in Ref. ^3] f° r 
the many-body inelastic piston problem. 

A particular solution to these macroscopic equations 
is symmetric cooling of both light particles, E\{t) — 

E 3 (t) = E Q [1 + V2#o(l - r 2 ) t]~ 2 , with E the ini- 
tial light particle energy, while the piston remains at 
Xiit) = 1/2. However, linear perturbation analysis shows 
that any small disturbance from symmetry grows and the 
piston is driven toward one of the walls, with an oscil- 
latory modulation that is periodic on logarithmic time 
scale |T^ . 

A typical piston trajectory that is obtained by numer- 
ically solving Eqs. ©-© with a slightly asymmetrical 
state is shown in Fig. ^ The numerical solution to the 
macroscopic equations and the simulation results for the 
three-particle system are extremely close over the time 
range 10 3 < t < 10 5 . However, after approximately 
10 5 time steps (for the case m 2 = 100 and r = 0.9), 
the macroscopic equations predict that inelastic collapse 
occurs, after which the piston sticks to one of the walls 
[l2|. In contrast, for the three-particle system, the piston 
localizes near one wall but continues to undergo small- 
amplitude, nearly regular oscillations on a logarithmic 
time scale (Fig.|2J. 

To help understand this discrepancy between the 
macroscopic approach and the simulations results for the 
three-particle system in the long-time limit, it is helpful 
to reconsider the elastic case r = 1. Here Eqs. iJJJJ and 
P)l can be immediately integrated, and substituting the 
results of these integrations into (JSJ gives 

d 2 x 2 _ Ai _ A 3 
7712 dt 2 ~ x\ (1 - X2) 3 ' 
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where are constants. This equation of motion de- 
scribes the oscillations of a particle in the effective poten- 
tial well V cS (x) = \{Aix~ 2 + A 3 (l - x)~ 2 }. This effective 
potential can be derived rigorously in the limit m 2 — » oo 
(see |l5| and also the Appendix of HR). Thus the long- 
time extreme excursions in the elastic system, which are 
not described by the effective potential, appear to stem 
from the fmiteness of the piston mass. 

By analogy, we anticipate that the macroscopic equa- 
tions JHJl-© should describe the final state for the inelas- 
tic piston in the m 2 — ► 00 limit. On the other hand, the 
log-periodic state emerges only when the piston mass is 
finite. This feature seems to play a parallel role as in the 
elastic system, in that departures from the predictions 
of the macroscopic equation arise only when the piston 
mass is finite. 



IV. EFFECTIVE RESTITUTION COEFFICIENT 

To understand the properties of the log-periodic oscil- 
lations, we map the three-particle system onto an equiva- 
lent two-particle system, from which the basic character- 
istics of the log-periodic state follow. The first step is to 
determine the net effect of the sequence of rattling colli- 
sions between the piston and a light particle as the piston 
approaches a wall and is ultimately reflected. We show 
in the Appendix that this collision sequence is equivalent 
to a f-body problem in which the piston is reflected from 
the wall with an effective restitution coefficient r c j(m2, r) 
that is smaller than the bare restitution coefficient r. 

Next, we exploit the symmetry breaking, in which the 
piston localizes near one wall, to reduce the initial three- 
body problem into an effective two-body problem that 
consists of the piston and one light particle. In this re- 
duced system, the piston collides inelastically with the 
wall with restitution coefficient r e g, while the light par- 
ticle collides elastically with the piston and inelastically 
with the other wall with restitution coefficient r (Fig. HJJ). 
Using this equivalence, we will determine the properties 
of the log-periodic state. 

For the initial step of determining the effective resti- 
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FIG. 6: Schematic space-time diagram of the particle trajec- 
tories in the log-periodic state (left) and the effective trajec- 
tories (right). 
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FIG. 7: (Color online) The effective restitution coefficient as 
a function of r for piston masses ni2 = 2, 5, fO, 30, fOO, 
and fOOO. The initial velocities are («i,«2) = (0, — f). The 
curves are the theoretical predictions from Eq. (lA.f 2L and 
the symbols correspond to simulation results. 



tution coefficient as a function of r and m 2 , the calcu- 
lational details are given in the Appendix and the final 
result for r e g(m 2 ,r) is quoted in Eq. (|A.f 2(1 . As shown 
in Fig.[7| r e fi decreases as r decreases and goes to zero as 
r approaches a critical value r c (771.2), quoted in Eq. (|A.7|I . 
that signals inelastic collapse. When r < r c , the effec- 
tive restitution coefficient is zero, and the result of the 
rattlings between the piston and the intervening light 
particle is inelastic collapse. For fixed r, notice also that 
r e ff decreases rapidly as m 2 is increased. 

With this effective restitution coefficient equivalence, 
we now reduce the original three-particle system to the 
equivalent two-particle system. Without loss of general- 
ity, we assume that the piston is close to the wall at x = 1 . 
The effective system then consists of a light particle at 
x± and the piston at X2, with < x\ < x 2 < 1 ■ For suffi- 
ciently large piston mass, the sequence of collisions con- 
sists of: (i) the piston making an effective collision with 
the wall with restitution coefficient r e g, (ii) the second 
light particle undergoing a bare inelastic collision with 
the other wall, and (hi) an elastic particle-piston collision, 
with each of these steps being non-overlapping. From the 
collision rules for each of these steps [see Eqs. QA.1J1 - 
(|A.2(l ]. the new velocities after each such collision se- 
quence are given in terms of the incoming velocities by 

/ m 2 - 1 2m 2 \ 

T T ff 1 

fv[\ I M M ( Vl \_ 

Ur 2 m2 -i W =Lv ' (6) 

V -M r -^T r ^l 

where M = 1 + m 2 , and vi,v[ < and v 2 , v' 2 > 0. The 
velocity vector after n such cycles is v™ = L"v. Diago- 
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nalizing L, we find (see also the Appendix) 

A^(A_-a)-A"(A+-a) 



,(«) 



A_ 

b(A1 - A") 

A_ - A+ 
(A" - A1)(A+ 



A, 



i'2 



a)(A_ 



6(A_-A + ) 
A I 1(A_ - a)- A"(A+ - a) 



"1 



A. 



A, 



-?' 2 



(7) 



(8) 



where a = r(m 2 — and 6 = — 2r c ffTO 2 /Af are the 

elements in the first row of L, and A± are the eigenvalues 
of matrix L, 



A, 



m 2 



1 



2M 



(r-r eS ) 1 ± 



Arr cS M 2 



(m 2 - l) 2 (r - r off ) 2 y ' 
(9) 

< 0, A i > 0, and 



Both eigenvalues are real, with A_ 
|A_| < |A+| < 1. 

We test this effective description of the collision dy- 
namics by comparing the exact piston trajectory in the 
three-particle system for a given m 2 and r with the pis- 
ton trajectory in the reduced two-particle system. After 
shifting the effective trajectory by an overall phase factor, 
both systems have visually indistinguishable periodic be- 
havior, thus confirming the validity of the coarse-grained 
approach. 



We now determine the relation that v[ n ^ and v^' 1 ' must 
satisfy for the effective two-particle system to be in a 
log-periodic state. For such a periodicity, successive col- 
lisions between particle 1 and the piston must occur at 
the same position xq for all n. Thus the time St^ for 
particle i to go from xq to its respective wall and return 
to xo in the n th cycle must be the same for both particles. 
That is, 



,(«) 



St[ n) = 



Thus 



Xq 1+T 



= 5ti n) = 



1 - Xq 1 + r e ff 



,(n). 



= 5t (7 



Tett 



\vj n) \ = l-tto r{\ + r cff ) 
|^")| xo r cS {l+r) 



(10) 



is a constant that is independent of n in the log-periodic 
state. Therefore \v ( 2 n+1) \/\v { ( l+1) \ = |« 2 " ) |/bi" ) |- Then 
using v( n+1 ) = Lv'™', we express v(" +1 ) in terms of v^ n -* 
and thereby obtain 



|^ n) | = (l + m 2 )A + - (mg - j> 
lu^l ~ 2m 2 r cff 



(11) 



where A + (m 2 , r) is the larger eigenvalue of L. Comparing 
Eqs. (jTUjl and l(TT|) finally yields 



Xq = 1 



(1 + r) [(1 + m 2 )A+ - (m 2 - l)r] 
2m 2 r(l + r off ) 



(12) 



and the amplitude of the piston oscillations is then 
A(rri2, r) = 1 — xq. This result agrees with the simulation 
results shown in Fig. even close to inelastic collapse. 
We may also compute the logarithmic period of the pis- 



ton oscillations. Since St^ n+1 ^/St^' 



o = i4")i 



we 



express |u 2 | in terms of \v[ n | and \v 2 | from Eq. JSJ, 
and then use Eq. I|ll|) to obtain 



~5tFT 



A, 



m 2 - 1 

1 

1 + m 2 



r e s [r - 



m 2 



1 



1 + m 2 



-A. 



> 1 



(13) 



Thus in the log-periodic state 8t^ n ' grows exponentially in 
the number of cycles n, as seen in our simulations. From 
the logarithmic period A(m 2 ,r) introduced in Sec. Ill Bl 
wehavei„_i = e~ A t n , so that St^ = t n (l-e~ A ). This 
relation then gives A = \n(t n+1 /t n ) = ln(<5^" +1 ) / 5t^) . 
The agreement between this prediction for A with simu- 
lation results (not shown) is again extremely good. 

Finally, the robustness of the log-periodic state can 
be understood in simple terms. Starting with an arbi- 
trary (not log-periodic) initial state, it is easy to show 
from Eqs. and (JSJ) that both v± and converge 
exponentially quickly with n to a state where the ratio 

\ v 2 \/\ v i I satisfies the condition (|ll|l that signals log- 
periodicity. This convergence occurs because for r c ff < r, 
(A_/A + ) n quickly goes to zero as n increases. In this 
sense, the log-periodic state is an attractor of the dy- 
namics. 



V. SUMMARY 

We investigated the dynamics of a three-particle sys- 
tem on the unit interval in which a massive particle (cor- 
responding to a piston) lies between 2 light particles. The 
particles collide elastically with the piston, but inelasti- 
cally with the walls. This toy model is meant to mimic 
the behavior of the inelastic piston problem in which 
a massive piston separates two inelastic gases, each of 
which contains many particles. The dynamics of this 
many-body problem is extremely rich. The piston moves 
nonmonotonically at early times and correspondingly the 
response of the two gases is also nonmonotonic. Even- 
tually there is an inelastic collapse in which one of the 
gases is compressed into a solid by the piston. 

One of the motivations for our study of the three- 
particle system was to capture some of the intriguing phe- 
nomenology of the many-particle inelastic piston prob- 
lem. A new feature of the three-particle system, how- 
ever, is that the piston settles into a log-periodic state at 
long times over a wide range of restitution coefficients, 
in which the period is constant on a logarithmic time 
scale. The characteristics of this log-periodic state can 
be understood in terms of a simple effective picture in 
which the rattling collision sequence between the piston 
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and the trapped light particle is replaced by an effec- 
tive inelastic collision between the piston and the wall, 
with effective restitution coefficient r Q g < r. This equiv- 
alence provides a satisfyingly complete account of the 
log-periodic state. Finally, it should be noted that the 
log-periodic behavior is a consequence of the finiteness 
of the piston mass. As to 2 increases, the amplitude of 
the oscillations decreases and as to 2 — > oo the inelastic 
collapse of the many-particle inelastic piston problem is 
recovered. 
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APPENDIX: CALCULATION OF THE 
EFFECTIVE RESTITUTION COEFFICIENT 



We use a matrix approach to compute the effective 
restitution coefficient r e g that describes the final veloc- 
ity of the piston at the end of the rattling collisions as 
a function of to.2 and the bare restitution coefficient r. 
Without loss of generality, we assume that the piston 
compresses particle 1 which then undergoes the rattling 
collision sequence. For concreteness, the light particle is 
taken to be at rest at x\ > with a wall at x — 0. A 
massive particle (the piston) approaches the light particle 
from the right with — — 1. Collisions between particle 
1 and the wall at x — are inelastic, with restitution 
coefficient r, while 1-2 collisions are elastic. After the 
rattling collision sequence ends, the piston recedes from 
the wall with velocity v' 2 = — «2 ^eff- 

The velocities after each collision are given in terms of 
the velocities before the collision by 



1 — m,2 
M 



M V1 



for the 1-2 collisions, and 



vi + 
m 2 



2m 2 
- 1 



V2 



M 



V2 , 



(A.l) 



-v\ r , 



V2 ■ 



(A.2) 



for the wall collision, where M = 1 + m 2 . Thus the com- 
bined effect of a 1-2 and an ensuing particle-wall collision 
is given by the composition of the two transformations 



implicit in Eqs. ljA.l|) and HA.2jl , Therefore 
1 — mo 



(- 



V 



M 
2 

M 



1 2 2m2 
r — ■ 



M 

7712 — 1 

M 



) 



JM eMv.(A.3) 



The velocity vector after n such cycles is given by 
v («) = mV°) , with v<°) = 




-1 



The collision sequence ends when the velocities of the 
two particles after n cycles satisfy — < 0, corre- 
sponding to the two particles receding from the wall with 
the piston moving faster than the light particle. We de- 
fine this situation as "escape" of the piston. The number 
of collisions for escape to occur is given by the small- 
est value of n that leads to the above conditions on the 
outgoing velocities. The effective restitution coefficient 
is then given by when n equals its value at escape. 

To determine the threshold value of n, we use the fact 
that (see, e.g., [T^ l 



(A.4) 



where S is the similarity matrix that diagonalizes M, 
and Mdiag = M n S is the diagonalized form of the 
transformation matrix. The eigenvalues of M are X± = 
(T±VT 2 - 4D)/2, where T = (m 2 - 1) (1 + r)/M is the 
trace and D — r is the determinant of M, 



(lW3 -l)(l + r) l_ 
2M I V 



4rM 2 



(m 2 - 1) 2 (1 + r) 2 J ' 
. ( A - 5) 

Consequently the similarity transformation matrix is 




where a — r(m 2 — ^)/M and b = —2iri2r/M are the el- 
ements of the first row of M, i.e., the matrix S consists 
of the eigenvectors of M arranged columnwise. Conse- 
quently S _1 = |Sp 1 S^, where |S| is the determinant of 
S, and is its transpose. 

Assembling these results, the velocity after n cycles 
(and 2n individual collisions) is 



» _ 



b X+ ' X ~ 



Xl(X+-a) A'l(A_- a ; 
V A_ - A i 



(n)\ 

.(«) 



(A.6) 



A_ - A, 



In the case where escape of the piston requires n + 1 
particle-particle collisions and n particle- wall collisions, 
we should multiply the transformation matrix M™ on the 
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left by the matrix defined by Eq. I|A.1|) to account for 

this last particle-particle collision. However, to compute 

only the final velocity of particle 2, it suffices to calculate 
v (n+i) from Eq 

Depending on the sign of the discriminant T 2 — AD, 
the eigenvalues A± can be real or complex. For r greater 
than a threshold value r c (m2), T 2 < AD. Thus A± are 
complex conjugates (note, however, that v^™) has always 
real components). At the threshold, T 2 = AD, leading 
to A + = A_, so that is undetermined. This inde- 
terminacy signals inelastic collapse: for r < r c (m 2 ) there 
is an infinite number of collisions in a finite time, and 

v[ n} ~ «4 U) > V "- The condition T 2 = AD gives the 
critical restitution coefficient for inelastic collapse: 



r c (m 2 ) 



(1 + ma)(l + m2- + Am 2 

(m 2 - l) 2 



(A.7) 



Notice that r c ~ 1 — 4/^/to2 in the limit of large mi- 
Equivalently, the condition T 2 — AD defines a critical 
mass Heir), such that inelastic collapse occurs for n%2 > 
Hc(r). We now find 

A«c(r)= (1 + r)(1 +:+^ F)+4r . (A.8) 

Note that for r close to 1, fi c (r) ~ 16(1 — r)~ 2 . 

For r > r c {m,2) the piston eventually escapes with ve- 
locity v 2 ™°' ) , where n is the number of cycles until escape. 

To determine no, define f(n) = u£ — i>2 • Initially 
/(0) = 1, and f(n) decreases as n increases and eventu- 
ally changes sign. Next, we define the real variable z by 



the condition f(z) — 0. From Eq. ljA.61 



/(*) = 



b{X z + - Af.) - A* (A+ -a) + X z _ (A_ - a) 

A_ - A, 



= 0. 
(A.9) 

Since A± are complex conjugates, we write X± — Q e & 
so that Eq. (|A.9|I becomes, using AfL — X z _ = 2iQ z sin(z/3) 
and a + b = —r, 



r sm(z/3) + Q sm[(z + l)/3] = , 



with solutions 



z(k) 



kn — tan 1 



1 sin j3 



Q cos/3 



(A.10) 



(A.ll) 



where k can be any integer number. The first solution 
that has a physical meaning (i.e., z > 1) corresponds 
to k = 1, so z = z(l). The number of collision cycles 
before escape is thus no = \z~\ , where \z~\ is the next 
integer larger than z. The escape velocity is v^f 1 "^ and 

r e s = «2 ° • However, for large enough piston mass, the 
number of collisions before escape is typically large, and 
we can approximate n by z. Hence, we finally obtain 
for the effective restitution coefficient, 



r e s(m 2 ,r) 



A^(A+-a)-Ai(A_ - a) 



A_ - A, 



(A.12) 



with z given by Eq. (|A.11|) with k = 1. A plot of r c s as 
a function of r is given in Fig. 
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